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Abstract. In this work we used Particle-In-Cell simulations to study the interaction of circularly polarised Alfven 
waves with one dimensional plasma density inhomogeneities transverse to the uniform magnetic field (phase 
mixing) in collisionless plasmas. In our preliminary work we reported discovery of a new electron acceleration 
mechanism, in which progressive distortion of the Alfven wave front, due to the differences in local Alfven speed, 
generates an oblique (nearly parallel to the magnetic field) electrostatic field. The latter accelerates electrons 
through the Landau resonance. Here we report a detailed study of this novel mechanism, including: (i) analysis 
of broadening of the ion distribution function due to the presence of Alfven waves and (ii) the generation of 
compressive perturbations due to both weak non-linearity and plasma density inhomogeneity. The amplitude 
decay law in the inhomogeneous regions, in the kinetic regime, is demonstrated to be the same as in the MHD 
approximation described by Heyvaerts and Priest (1983). 
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1. Introduction 

The study of the interaction of Alfven waves (AWs) with 
plasma inhomogeneities is important for both astrophys- 
ical and laboratory plasmas. This is because both AWs 
and inhomogeneities often coexist in these physical sys- 
tems. AWs are believed to be good candidates for plasma 
heating, energy and momentum transport. On the one 
hand, in many physical situations AWs are easily ex- 
citable (e.g. through convective motion of the solar in- 
terior) and thus they are present in a number of astro- 
physical systems. On the other hand, these waves dissi- 
pate due to the shear viscosity as opposed to compres- 
sive fast and slow magnetosonic waves which dissipate due 
to the bulk viscosity. In astrophysical plasmas shear vis- 
cosity is extremely small as compared to bulk viscosity. 
Hence, AWs are notoriously difficult to dissipate. One of 
the possibilities to improve AW dissipation is to introduce 
progressively decreasing spatial scales, SI 0, into the 
system (recall that the classical dissipation is cx 6l~-^). 
Heyvaerts and Priest have proposed (in the astrophysi- 
cal context) one such mechanism, called AW phase mix- 



ing I Hevvaerts fc Priestlll983|) . It occurs when a linearly 
polarised AW propagates in the plasma with a one di- 
mensional density inhomogeneity transverse to the uni- 
form magnetic field. In such a situation the initially plane 
AW front is progressively distorted because of different 
Alfven speeds across the field. This creates progressively 
stronger gradients across the field (in the inhomogeneous 
regions the transverse scale collapses to zero), and thus 
in the case of finite resistivity, dissipation is greatly en- 
hanced. Hence, it is believed that phase mixing can pro- 
vide significant plasma heating. Phase mix ing could be 
also im portant for laboratory plasmas. .Haseeawa fc: CherJ 
proposed the heating of collisionless plasma by 
utilising spatial phase mixing by shear Alfven wave res- 
onance and discussed potential applications to toroidal 
plasma. A significant amount of work has been done in 
the context of heating open magnetic structu res in the so- 
lar corona llHevva,erts fc p'riest. 1983 Nocera et alJIlflS 
Parked ll99U iNakariakov et all ll997L__lDeMoortel et a 
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2000, Botha et al 2000, Tsiklauri et al 2001, Hood et al.' 
2002, Tsiklauri fc Nakariako v 2002, Tsiklauri ct al. 2002, 
2003). All phase mixing studies so far have been performed 
in the MHD approximation, however, since the transverse 
scales in the AW collapse progressively to zero, the MHD 
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approximation is inevitably violated. This happens when 
the transverse scale approaches the ion gyro-radius 
and then electron gyro-radius r,,. Thus, we proposed to 
study the phase mixing effect in the kinetic regime, i.e. 
we go beyond a MHD approximation. P reliminary results 
were reported in lTsiklauri et all l)2005|) . where we discov- 
ered a new mechanism for the acceleration of electrons 
due to wave-particle interactions. This has important im- 
plications for various space and laboratory plasmas, e.g. 
the coronal heating problem and acceleration of the so- 
lar wind. In this paper we present a full analysis of the 
discovered effect including an analysis of the broadening 
of the ion distribution function due to the presence of 
Alfven waves and the generation of compressive perturba- 
tions due to both weak non-linearity and plasma density 
inhomogeneity. 



2. The model 

We used 2D3V, the fully relativistic, electromagnetic, 
particle-in-cell (PIC) code with MPI para llelisation, 
modified f rom the 3D3V TRISTAN code l|Bunemanl 
11993. D.67|) . The system size is = 5000 A and Ly = 
200A, where A(= 1.0) is the grid size. The periodic 
boundary conditions for x- and y-directions are imposed 
on particles and fields. There are about 478 million elec- 
trons and ions in the simulation. The average number of 
particles per cell is 100 in low density regions (see be- 
low). The thermal velocity of electrons is Vth,e = 0.1c 
and for ions is Vth i — 0.025c. The ion to electron mass 



ratio is rrii/me = 16. The time step is u>peAt = 0.05. 
Here iVpe is the electron plasma frequency. The Debye 
length is vth.e/^pe = 1-0. The electron skin depth is 
c/u>pe = lOA, while the ion skin depth is c/upi = 40A. 
Here LOpi is the ion plasma frequency. The electron Larmor 
radius is Vth,e/^ce — 1-OA, while the same for ions is 
Vth,i/^ci — 4. OA. The external uniform magnetic field, 
Bq(= 1.25), is in the x-direction and the initial electric 
field is zero. The ratio of electron cyclotron frequency to 
the electron plasma frequency is ojce/ijJpe = 1-0, while the 
same for ions is cod/uipi = 0.25. The latter ratio is essen- 
tially Va/c~ the Alfven speed normahsed to the speed of 
light. Plasma (3 = 2{u}pe/uJceY{vth,elcY = 0.02. Here ah 
plasma parameters are quoted far away from the density 
inhomogeneity region. The dimensionless (normalised to 
some reference constant value of no = 100 particles per 
cell) ion and electron density inhomogeneity is described 

by 



ni{y) = ne{y) = 1 -|- 3exp 



y- IOOAV 



50A 



F{y). (1) 



This means that in the central region (across the y- 
direction), the density is smoothly enhanced by a factor 
of 4, and there are the strongest density gradients hav- 
ing a width of about 50 A around the points y = 5 1.5 A 
and y = 148. 5A. The background temperature of ions and 



electrons, and their thermal velocities are varied accord- 
ingly 



T,iy)/To^T,{y)/To^Fiy)-\ 

Vth,t/VtO = Vth.e/VeO = ^(j/)~^^^, 



(2) 
(3) 



such that the thermal pressure remains constant. Since the 
background magnetic field along the x-coordinate is also 
constant, the total pressure remains constant too. Then 
we impose a current of the following form 



dtE, 



dtE, 



-Jo sm(ujdt) (1 - exp [-{t/to)^] ) 



'jocos(wdt) (l 



exp 



(4) 
(5) 



Here ujd is the driving frequency which was fixed at 
Ud = O.Swci- This ensures that no significant ion-cyclotron 
damping is present. Also, dt denotes the time derivative, to 
is the onset time of the driver, which was fixed at bO/ujpe 
i.e. 3.125/wci. This means that the driver onset time is 
about 3 ion-cyclotron periods. Imposing such a current 
on the system results in the generation of left circularly 
polarised AW, which is driven at the left boundary of sim- 
ulation box and has spatial width of 1 A. The initial ampli- 
tude of the current is such that the relative AW amplitude 
is about 5 % of the background (in the low density homo- 
geneous regions), thus the simulation is weakly non- linear. 

3. Main results 

Because no initial (perpendicular to the external mag- 
netic field) velocity excitation wa s imposed in addition to 
the ab o ve specified curr e nts (c f. 
(|2002|) iDelZanna et alJ l)200: 
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the circularly polarised AW excited (driven 
at the left boundary is split into two circularly polarised 
AWs which travel in opposite directions. The dynamics of 
these waves as well as othe r physical quant i ties is shown 
in Fig. 1. (cf. Fig. 1 from iTsiklauri et "all l)2005(l where 
Bz and Ey, the circularly polarised Alfven wave compo- 
nents, were shown for three different times). A typical 
simulation, untill t — 875/ujce = 54.69/tJci takes about 
8 days on the parallel 32 dual 2.4 GHz Xeon processors. 
It can be seen from the figure that because of the peri- 
odic boundary conditions, a circularly polarised AW that 
was travelling to the left, reappeared on the right side of 
the simulation box. The dynamics of the AW {Bz,Ey) 
progresses in a similar manner as in the MHD, i.e. it 
phase mixes l|Hevvaerts fc Pries il ll983l) . In other words, 
the middle region (in y-coordinate) , i.e. 50A < y < 150A, 
travels slower because of the density enhancement (note 
that VA{y) oc I/a/ ni{y)). This obviously causes a dis- 
tortion of initially plain wave front and the creation of 
strong gradients in the regions around y « 50 and 150. 
In the MHD approximation when resistivity, 77, is finite, 
the AW is strongly dissipated in these regions. This effec- 
tively means that the outer and inner parts of the trav- 
elling AW are detached from each other and propagate 
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independently. This is why the effect is called phase mix- 
ing - after a long time (in the case of developed phase 
mixing), phases i n the wave front becom e effectively un- 
correlated. Before ITsiklauri et alJ l|2005|) . it was not clear 
what to expect from our PIC simulation. The code is 
coUisionless and there are no sources of dissipation in it 
(apart from the possibility of wave-particle interactions). 
It is evident from Fig. 1 that in the developed stage of 
phase mixing {t = 54.69 /ujd), the AW front is substan- 
tially damped in the strongest density gradient regions. 
Contrary to the AW {B^ and Ey) dynamics we do not see 
any phase mixing for By and E^. The latter two behave 
similarly. It should be noted that E^ contains both driven 
(see Eq.(4) above) and non-linearly generated fast mag- 
netosonic wave components, with the former being domi- 
nant over the latter. Since B^ is not driven, and initially 
its perturbations are absent, we see only non-linearly gen- 
erated slow magnetosonic perturbations confined to the 
regions of strongest density gradients (around y 50 and 
150). Note that these also have rapidly decaying ampli- 
tude. Also, we gather from Fig. 1 that the density per- 
turbation (« 10%) which is also generated through the 
weak non-linearity is present too. These are propagating 
density oscillations with variation both in overall magni- 
tude (perpendicular to the figure plane) and across the y- 
coordinate, and they are mainly confined to the strongest 
density gradients regions (around y « 50 and 150). Note 
that dynamics of B^ , By and B^ with appropriate geo- 
metrical switching (because in our geometry the uniform 
magnetic field lies a long the x-co ordinate) is in qualita- 
tive agreement with iBotha et all l|2000 'l (cf their Fig. 9). 
The dynamics of remaining E^ component is treated sep- 
arately in the next figure. It is the inhomogeneity of the 
medium ni{y), VA{y), i.e. d/dy ^ 0, is the cause of weakly 
non- linear coupling of t he AWs to the compres s ive ni odes 
(see INakariakov et alJ l(l997l) . iTsikla^iri et aD ((200 ll) for 
further details). 

In Fig. 2 we try to address the question of where the 
AW energy went? (as we saw strong decay of AWs in the 
regions of strong density gradients). Thus in Fig. 2 we 
plot Ex, the longitudinal electrostatic field, and electron 
phase space {Vx/c vs. x and Vx/c vs. y) for different times. 
In the regions around y « 50 and 150, for later times, 
a significant electrostatic field is generated. This is the 
consequence of stretching of the AW front in those re- 
gions because of the difference in local Alfven speed. In 
the middle column of this figure we see that exactly in 
those regions where Ex is generated, many electrons are 
accelerated along x-axis. We also gather from the right 
column that for later times (t = 54:.69/LUci), the number 
of high velocity electrons is increased around the strongest 
density gradient regions (around y « 50 and 150). Thus, 
the generated Ex field is somewhat oblique (not exactly 
parallel to the external magnetic field). Hence, we con- 
clude that the energy of the phase-mixed AW goes into 
acceleration of electrons. Line plots of Ex show that this 
electrostatic field is strongly damped, i.e. the energy is 
channelled to electrons via Landau damping. 



In Fig. 3 we investigate ion phase space {Vz/c vs. x and 
Vz/c vs. y) for the different times. The reason for choice of 
Vz will become clear below. We gather from this plot that 
in Vz/c vs. X phase space, clear propagating oscillations 
are present (left column). These oscillations are of the in- 
compressible, Alfvenic "kink" type, i.e. for those xs where 
there is an increase of greater positive velocity ions, there 
is also a corresponding decrease of lower negative velocity 
ions. In the Vz/c vs. y plot we also see no clear acceleration 
of ions. 

We next look at the distribution functions of electrons 
and ions before and after the phase mixing took place. 
In Fig. 4 we plot distribution functions of electrons and 
ions at t = and t = 54:.69/uJci- Note that at < = the 
distribution functions do not look as purely Maxwellian 
because of the fact that the temperature varies across 
the y-coordinate (to keep total pressure constant) and 
the graphs are produced for the entire simulation domain. 
Also, note that for electrons in /(T4) there is a substantial 
difference at t — 54.69/cJci to its original form because of 
the aforementioned electron acceleration. We see that the 
number of electrons having velocities 14 = ±(0.1 — 0.3)c 
is increased. Note that the acceleration of electrons takes 
place mostly along the external magnetic field (along the 
x-coordinate). Thus, very little electron acceleration oc- 
curs for Vy or Vz (solid and dotted curves practically over- 
lap each other). For the ions the situation is different: we 
see broadening of the ion velocity distribution functions 
in Vz and Vy (that is why we have chosen to present the 
Vz component of ion phase space in Fig. 3). 

The reason for this broadening of the ion distribution 
function becomes clear in Fig. 5 where we plot kinetic 
energy x, y, z components (oc Vxy z) ^^^d total kinetic en- 
ergies for electrons (top row) and ions (bottom row). For 
ions we gather that y and z components of the kinetic en- 
ergy (bottom left figure) oscillate in anti-phase and their 
oscillatory part perfectly cancels out in the total energy 
(bottom right figure). Thus, the broadening of the y and 
z components of the ion velocity distribution functions 
is due to the presence of AWs (usual wave broadening, 
which is ac tually observed e . g. in t he solar corona an d 
solar wind, iBe lcher fc Davij ((l97ll) . ISmith et all l(l995l) . 
iBavassano et al. (.2000) )■ and hence there is no ion accel- 
eration present. Note that y and z components and hence 
total kinetic energy of ions is monotonously increasing due 
to continuous AW driving. Note that no significant motion 
of ions along the field is present. For electrons, on the other 
hand, we see a significant increase of the x component 
(along the magnetic field) of kinetic energy which is due 
to the new electron acceleration mechanism discovered by 
us (cf. Fig. 10). Note that for ions the y component reaches 
lower values (than the z component) because of lower AW 
velocity in the middle part of the simulation domain. 

The next step is to check whether the increase in elec- 
tron velocities really comes from the resonant wave parti- 
cle interactions. For this purpose in Fig. 6, left panel, we 
plot two snapshots of the Alfven wave Bz{x,y — 148) 
component at instances t = 54:.69/ujci (solid line) and 
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Fig. 1. Contour (intensity) plots of electromagnetic field components and electron density at time t — 54:.Q9/uJci 
(developed stage of phase mixing). The phase mixed Alfven wave components are and Ey. The excitation source 
is at the left boundary. Because of periodic boundary conditions, the left-propagating AW re-appears from the right 
side of the simulation box. Note how the (initially plain) AW is stretched because of differences in local Alfven speed 
across the y-coordinate. Significant (w 10%) density fluctuations can be seen. 



t = 46.87 /ojci (dotted line). The distance between the two 
upper leftmost peaks (which is the distance travelled by 
the wave in the time span between the snapshots) is about 
6L = 150A = 15{c/ujpe). The time difference between the 
snapshots is St = 7.82/uJci- Thus, the measured AW speed 
at the point of the strongest density gradient {y — 148) 
is = SL/St = 0.12c. We can also work out the Alfven 
speed from theory. In the homogeneous low density re- 
gion the Alfven speed was set to be Va{oo) — 0.25c. From 
Eq.(l) it follows that for y = 148 the density is increased 
by a factor of 2.37 which means that the Alfven wave speed 
at this position is Va(148) = 0.25/a/2^c = 0.16c. The 
measured (0.12c) and calculated (0.16c) Alfven speeds in 
the inhomogeneous regions do not coincide. This is prob- 
ably because the AW front is decelerated (due to momen- 
tum conservation) as it passes on energy and momentum 
to the electrons in the inhomogeneous regions (where elec- 
tron acceleration takes place). However, this possibly is 
not the case if wave-particle interactions play the sa me role 
as dissipation in MHD (Saedeev fc GalceA^ll969|) : Then 



wave-particle interactions would result only in the de- 
crease of the AW amplitude (dissipation) and not in its 
deceleration. If we compare these values to Fig. 4 (top left 
panel for f{Vx)), we deduce that these arc the velocities 
> 0.12c above which electron numbers with higher veloc- 
ities are greatly increased. This deviation peaks at about 
0.25c which, in fact, corresponds to the Alfven speed in 
the lower density regions. This can be explained by the 
fact the electron acceleration takes place in wide regions 
(cf. Fig. 2) along and around y « 148 (and y ~ 51) - 
hence the spread in the accelerated velocities. In Fig. 6 
we also plot a visual fit curve (dashed line) to quantify 
the amplitude decay law for the AW (at t = 54.69/a;ci) 
in the strongest density inhomogeneity region. The fitted 

(dashed) cure is represented by 0.056 exp — (a;/1250)^ . 
There is a surprising si milarity of this fit w i th th e MHD 
approximation results. iHevvaerts fc Priest l)l983|) found 
that for large times (developed phase mixing), in the case 
of a harmonic driver, the amplitude decay law is given 
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Fig. 2. Left column: contour plots of the generated electrostatic fielcl nearly parallel to the external magnetic 
field at instances: t = (0, 31.25, 54. 69)/wci- Central column: T4/c versus x of electron phase space at the same times. 
To reduce figure size, only electrons with Vx > 0.15c were plotted. Right column: Vx/c versus y of electron phase 
space at the same times. Only electrons with Vx > 0.12c were plotted (note the dip in the middle due to the density 
inhomogeneity across ?/-coordinate) . 




Fig. 3. Left column: V^/c versus x of ion phase space at instances: t = (0, 31.25, 54.69) /uid- Right column: V^/c versus 
y of ion phase space at the same times. Only ions with Vz > 0.03c are plotted (note the dip in the middle due to the 
density inhomogeneity across the y-coordinate) . 



by (X exp — ^ 



6Vi 



which is much faster than the 

usual resistivity dissipation oc exp{—r]x). Here is the 
derivative of the Alfven speed with respect to the y- 
coordinate. The most interesting fact is that even in the 



kinetic approximation the same oc exp(— Aa;^) law holds 
as in MHD. 

In MHD, finite resistivity and Alfven speed non- 
uniformity are responsible for the enhanced dissipation via 
phase mixing. In our PIC simulations (kinetic phase mix- 
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Fig. 6. Left: two snapshots of the Alfven wave Bz{x,y = 148) component at instances t = 54.69 /ojd (solid line) 

and t — 46.87 /u!ci (dotted line). The dashed line represents fit 0.056 exp -(a;/1250) . Center: Bz{x,y = 10) (low 

density homogeneous region), Bz{x,y — 100) (high density homogeneous region). Note the differences in amplitudes 
and propagation speeds, which are consistent with the equilibrium density and thus Alfven speed dependence on 
2/-coordinate. 



ing), however, we do not have dissipation and collisions 
(dissipation). Thus, in our case, wave-particle interactions 
play th e same role as resistivit y rj in the MHD phase 
mixing ijSaedeev fc Galeevlll969|) . No significant AW dis- 
sipation was found away from the density inhomogene- 
ity regions (Fig. 6 middle and right panels, note also the 



differences in amplitudes and propagation speeds, which 
are consistent with the imposed density and hence Alfven 
speed variation across the y-coordinate). This has the 
same explanation as in the case of MHD - it is in the 
regions of density of inhomogeneities (V^ 7^ 0) that the 
dissipation is greatly enhanced, while in the regions where 
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Fig. 5. Top row: kinetic energies (calculated by x (solid), 
y (dotted), z (dashed) velocity (squared) components) 
of electrons (left), and total kinetic energy for electrons 
(right) as a function of time. Bottom row: As above but 
for ions. The units on y-axis are arbitrary. 



= there is no substantial dissipation (apart from the 
classical oc exp{—r]x) one). In the MHD approximation, 
the aforementioned amplitude decay law is derived from 
the diffusion equation, to which MHD e quations reduce 
for la rge times (developed phase mixing ijTsiklauri et alJ 
It seems that the kinetic description leads to the 
same type of diffusion equation. It is unclear at this stage, 
however, what physical quantity would play the role of re- 
sistivity rj (from the MHD approximation) in the kinetic 
regime. 

3.1. Homogeneous plasma case 

In order to clarify the broadening of the ion velocity dis- 
tribution function and also for a consistency check we 
performed an additional simulation in the case of ho- 
mogeneous plasma. Now the density was fixed at 100 
ions/electrons per cell in the entire simulation domain 
and hence plasma temperature and thermal velocities were 
fixed too. In such a set up no phase mixing should take 
place as the AW speed is uniform. 

In Fig. 7 we plot the only non-zero (above noise level) 
components at t — 54.69/ct;ci, which are left circularly po- 
larised AW fields: B^, By, E^, Ey. Note there are no B^ , E^ 
or density fluctuations present in this case (cf. Fig. 1) as it 
is the plasma inhomogeneity that facilitates the coupling 
between AW and the compressive modes. 

In Fig. 8 we plot ion phase space {Vz/c vs. x and Vz/c 
vs. y) in the homogeneous plasma case for different times 
(cf. Fig. 3). We gather from the graph that propagating, 
incompressible, Alfvenic "kink" type oscillations are still 
present (left column), while no significant ion accelera- 
tion takes place (right column) . This is better understood 
from Fig. 9 where we plot electron and ion distribution 
functions for t = Q and t = 54.69/cL'ci (as in Fig. 4) for 
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Fig. 10. As in Fig. 5 but for the case of homogeneous 
plasma density (no phase mixing). 

the homogeneous plasma case. There are three notewor- 
thy points: (i) no electron acceleration takes place because 
of the absence of phase mixing; (ii) there is (as in the inho- 
mogeneous case) broadening of the ion velocity distribu- 
tion functions (in Vy and Vz) due to the present AW (wave 
broadening); (iii) The distribution now looks Maxwellian 
(cf. Fig. 4) because the distribution function holds for the 
entire homogeneous region. 

In Fig. 10 we plot the kinetic energy x, y, z components 
(cx V^^^y^z) and total kinetic energy for electrons (top row) 
and ions (bottom row). For ions we see (as in the inho- 
mogeneous case) that y and z components of the kinetic 
energy (bottom left figure) oscillate in anti-phase and their 
oscillatory part perfectly cancels out in the total energy 
(bottom right figure). Thus, the broadening of the y and 
z components of the ion velocity distribution functions is 
due to the presence of AWs, and, in turn, there is no ion 
acceleration. Again the y and z components and hence the 
total kinetic energy of the ions is monotonously increas- 
ing due to continuous AW driving. No significant motion 
of ions along the field is present. For electrons we do not 
observe any acceleration due to the absence of phase mix- 
ing (cf. Fig. 5). Note that the y component now attains 
the same values as the z component (bottom left figure) 
because of the same AW velocity in the entire simulation 
domain. 

4. Discussion 

In our preliminary work l)Tsiklauri et alJ I200,'t() we out- 
lined the main results of our newly-discovered mechanism 
of electron acceleration. Here we presented a more detailed 
analysis of the phenomenon. We have established the fol- 
lowing: 

— Progressive distortion of the Alfven wave front, due to 
the differences in local Alfven speed, generates oblique 
(nearly parallel to the magnetic field) electrostatic 
fields, which accelerate electrons. 
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Fig. 7. As in Fig. 1 but for the case of homogeneous plasma density (no phase mixing) . Note only non-zero (above 
noise level) components are plotted. There are no B^, or density fluctuations present in this case. 
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Fig. 8. As in Fig. 3 but for the case of homogeneous plasma density (no phase mixing). 



— The amplitude decay law in the inhomogeneous re- 
gions, in the kinetic regime, is shown to be the 
same as in the MHD a pproximation described by 
iHevvaerts fc Priest! l)l983|) . 

— The density perturbations (« 10% of background) 
are generated due to both the weak non-linearity and 
plasma inhomogeneity. These are propagating density 
oscillations with variations both in overall magnitude 
and across the y-coordinate. They are mainly con- 
fined to the strongest density gradients regions (around 



J/ « 50 and 150) i.e. edges of the density structure (e.g. 
boundary of a solar coronal loop) . Longitudinal to the 
external magnetic field, B^, perturbations are also gen- 
erated in the same manner, but with smaller (« 3%) 
amplitudes. 

— Both in the homogeneous and inhomogeneous cases the 
presence of AWs causes broadening of the perpendicu- 
lar (to the external magnetic field) ion velocity distri- 
bution functions, while no ion acceleration is observed. 
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Fig. 9. As in Fig. 4 but for the case of homogeneous plasma density (no phase mixing). 



In the MED app roximation iHood et al.l (j2QQ2l) and 
iTsiklauri et al ] (|2003') showed that in the case of locahsed 
Alfven pulses, Heyvaerts and Priest's amplitude decay for- 
mula cx exp(— Aa;'^) (which is true for harmonic AWs) is 
replaced by the power law Bz oc x"^/^. A natural next 
step forward would be to check whether in the case of 
localised Alfven pulses the same power law holds in the 
kinetic regime. 

After this study was complete w e became aware of a 
study bv lVasauez fc Hollwed l|2004^ ■ who used a hybrid 
code (electrons treated as a neutralising fluid, with ion 
kinetics retained) as opposed to our (fully kinetic) PIC 
code, to simulate resonant absorption. They found that a 
planar (body) Alfven wave propagating at less than 90° to 
a background gradient has field lines which lose wave en- 
ergy to another s e t of fi eld lines by cross- field transport. 
Further, IVasaued l|2004 found that when perpendicular 
scales of the order of 10 proton incrtial lengths {IQc/ujpi) 
develop from wave refraction in the vicinity of the resonant 
field lines, a non-propagating density fluctuation begins to 
grow to large amplitudes. This saturates by exciting highly 
oblique, compressive and low-frequency waves which dis- 
sipate and heat protons. These processes lead to a faster 
development of small scales across the magnetic field, i.e. 
the phase mixing mechanism, studied here. 
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